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A new type of Ei Nino event has been identified in the last decade. During “warm pool" El 
Nino (WPEN) events, sea surface temperatures (SSTs) in the central equatorial Pacific are 
warmer than average. The El Nino signal propagates poleward and upward as large-scale 
atmospheric waves, causing unusual weather patterns and warming the polar stratosphere. 
In austral summer, observations show that the Antarctic lower stratosphere is several degrees 
(K) warmer during WPEN events than during the neutral phase of El Nino/Southern Oscillation 
(ENSO). Furthermore, the stratospheric response to WPEN events depends of the direction of 
tropical stratospheric winds: the Antarctic warming is largest when WPEN events are 

coincident with westward winds in the tropical lower and middle stratosphere i.e., the 
westward phase of the quasi-biennial oscillation (QBO). Westward winds are associated with 
enhanced convection in the subtropics, and with increased poleward wave activity. 

In this paper, a new formulation of the Goddard Earth Observing System Chemistry-Climate 
Model, Version 2 (GEOS V2 CCM) is used to substantiate the observed stratospheric response 
to WPEN events. One simulation is driven by SSTs typical of a WPEN event, while another 
simulation is driven by ENSO neutral SSTs; both represent a present-day climate. Differences 
between the two simulations can be directly attributed to the anomalous WPEN SSTs. During 
WPEN events, relative to ENSO neutral, the model simulates the observed increase in 
poleward planetary wave activity in the South Pacific during austral spring, as well as the 
relative warming of the Antarctic lower stratosphere in austral summer. However, the 
modeled response to WPEN does not depend on the phase of the QBO. The modeled 
tropical wind oscillation does not extend far enough into the lower stratosphere and upper 
troposphere, likely explaining the model’s insensitivity to the phase of the QBO during WPEN 
events. 
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1 


Abstract 


2 A new formulation of the Goddard Earth Observing System Chemistry-Climate Model, 

3 Version 2 (GEOS V2 CCM), with an improved general circulation model and an internally 

4 generated quasi-biennial oscillation (QBO), is used to investigate the response of the 

5 Antarctic stratosphere to (1) warm pool El Nino (WPEN) events and (2) the sensitivity of this 

6 response to the phase of the QBO. Two 50-year time-slice simulations are forced by 

7 repeating annual cycles of sea surface temperatures and sea ice concentrations composited 

8 from observed WPEN and neutral ENSO (ENSON) events. In these simulations, greenhouse 

9 gas and ozone-depleting substance concentrations represent the present-day climate. The 

10 modelled responses to WTEN, and to the phase of the QBO during WPEN, are compared with 

1 1 NASA’s Modern Era Retrospective-Analysis for Research and Applications (MERRA) 

12 reanalysis. 

13 

14 WPEN events enhance poleward planetary wave activity in the central South Pacific region 

15 during austral spring, leading to relative warming of the Antarctic lower stratosphere in 

16 November,' December. During the easterly phase of the QBO (QBO-E), the GEOS V2 CCM 

17 reproduces the observed 4-5 K warming of the polar region at 50 hPa, in the WPEN 

18 simulation relative to ENSON. 

19 

20 In the recent past, the response to WPEN events was sensitive to the phase of the QBO: the 

21 enhancement in planetary wave driving and the lower stratospheric wanning signal were 

22 mainly associated with WPEN events coincident with QBO-E. In the GEOS V2 CCM, 

23 however, the Antarctic response to WPEN events is insensitive to the phase of the QBO: the 

24 modelled response is always easterly QBO-like. OLR, streamfunction and Rossby wave 

25 energy diagnostics are used to show that the modelled QBO does not extend far enough into 

26 the lower stratosphere and upper troposphere to modulate convection and thus planetary wave 

27 activity in the south central Pacific. 
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1 1 Introduction 


2 Recent literature has identified two types of El Nino events. Conventional or “cold tongue” 

3 El Nino (CTEN) events are characterized by positive sea surface temperature (SST) 

4 anomalies in the eastern equatorial Pacific (Nino 3 region) (Rasmusson and Caipenter, 1982; 

5 Kug et al., 2009), while “warm pool” El Nino (WPEN) events are characterized by positive 

6 SST anomalies in the central equatorial Pacific (Nino 4 region) (Larkin and Harrison, 2005; 

7 Ashok et ah, 2007; Kug et ah, 2009). 

8 

9 The Northern Hemisphere (NH) stratosphere response to CTEN events is well recognized. El 

10 Nino-related warming of the Arctic stratosphere has been identified in both observational 

11 (Bronnimann et ah, 2004; Free and Seidel, 2009) and modelling studies (Sassi et ah, 2004; 

12 Manzini et ah, 2006; Cagnazzo et ah, 2009). This warming is a response to increased 

13 planetary wave driving: Garfinkel and Hartmann (2008) showed that the extra-tropical 

14 tropospheric teleconnections produced during El Nino events weaken the Arctic vortex, 

15 leading to higher stratospheric temperatures during the NH winter season. Furthermore, the 

16 phase of the quasi-biennial oscillation (QBO) modulates the stratospheric response to CTEN: 

1 7 The Arctic vortex is weakest in years when CTEN events coincide with the easterly phase of 

18 the QBO (Garfinkel and Hartmann, 2007). 

19 

20 Hurwitz et ah (2011) identified a robust response to WPEN in the Southern Hemisphere (SH) 

2 1 spring and summer. By compositing observed WPEN and ENSO neutral (ENSON) events, 

22 these authors used meteorological reanalyses to show the poleward extension and increased 

23 strength of the South Pacific Convergence Zone SPCZ during WPEN events, as compared 

24 with ENSON, in austral spring. This configuration of convective activity in the south Pacific 

25 favoured an enhancement of planetary wave activity in the upper troposphere, and on average, 

26 higher polar stratospheric temperatures and a weakening of the Antarctic polar jet in austral 

27 summer. Analogously to the NH response to CTEN events, the authors found the SH 

28 response to WPEN events to be dependent on the phase of the quasi-biennial oscillation 

29 (QBO): planetary wave driving was strongest in WPEN years coincident with the easterly 

30 phase of the QBO. 

31 
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1 While the conclusions reached by Hurwitz et al. (201 1) were based on physical arguments, 

2 statistically robust, and consistent amongst several reanalysis datasets, the scope of the 

3 authors’ study was limited by the small number of WPEN events during the satellite era. In 

4 particular, the WPEN and easterly QBO composite consisted of just three events. Hurwitz et 

5 al. (2011) concluded that simulations forced by repeating WPEN and ENSON boundary 

6 conditions would greatly increase the sample size and strengthen their findings. Such 

7 simulations require a model with a well-resolved stratosphere, a QBO, and because the 

8 stratospheric response to WPEN occurs in austral spring and summer, interactive polar ozone 

9 chemistry. 

10 

1 1 This paper presents the first attempt to simulate the SH stratospheric response to WPEN 

12 events. Section 2 describes the atmospheric datasets and time-slice simulations with the 

13 Goddard Earth Observing System Chemistry-Climate Model, Version 2 (GEOS V2 CCM) to 

14 be analyzed. In Section 3, OLR, tropospheric stationaiy wave patterns, a Rossby wave source 

15 diagnostic, eddy heat flux at 100 hPa and stratospheric temperatures illustrate the observed 

16 and modelled responses to WPEN events, as well as modulation of the WPEN response by the 

17 QBO. Section 4 provides a summary of the results and a brief discussion. 

18 

19 2 Model and data sources 

20 2.1 Atmospheric datasets 

21 The Modem Era Retrospective- Analysis for Research and Applications (MERRA) is used to 

22 calculate streamfunction, eddy heat flux, temperature, wind and planetary wave energy 

23 diagnostics for the 1979-2009 period. MERRA is a reanalysis dataset based on an extensive 

24 set of satellite observations and on the Goddard Earth Observing System Data Analysis 

25 System, Version 5 (GEOS-5)(Bosilovich et al., 2008; Rienecker et al., in preparation). The 

26 MERRA reanalysis has vertical coverage up to 0. 1 hPa, and for this study, is interpolated to 

27 1.25° x 1 .25° horizontal resolution. Hurwitz et al. (2011) showed that the response to WPEN 

28 was consistent amongst several meteorological reanalyses, including MERRA. 

29 

30 Zonal winds from the MERRA reanalysis are used to characterize the phase of the QBO. For 

31 each year between 1979 and 2009, November/December mean zonal winds at 50 hPa, 
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1 between 10°S and 10°N, form a QBO index: QBO easterly years (QBO-E) are identified 

2 when the QBO index is less than -2 m s' 1 ; QBO westerly years (QBO-W) are identified when 

3 the QBO index is larger than 2 m s' 1 . 

4 

5 Outgoing longwave radiation (OLR) between 1979 and 2009 is obtained from the NOAA 

6 interpolated OLR dataset (Liebmann and Smith, 1996), provided by NOAA/OAR/ESRL 

7 Physical Sciences Division, Boulder, CO. 

8 

9 WPEN and ENSON events are identified using SON seasonal mean SST anomalies in the 

10 Nino 3 and Nino 4 regions (www.cpc.noaa.gov/data/indices), as in Hurwitz et al. (2011). 

1 1 Because of the anomalous conditions in the SH stratosphere during the 2002 winter season 

12 (i.e., Newman and Nash, 2005), the 2002-2003 WPEN event has been excluded from the 

13 present analysis. The observed WPEN and ENSON events considered in the following 

14 analysis are listed in Table 1. Three WPEN events are coincident with QBO-E (WPEN/QBO- 

15 E) and three are coincident with QBO-W (WPEN/QBO-W). 12 ENSON events are 

1 6 distributed throughout the satellite era. 

17 

18 2.2 A new formulation of the GEOS V2 CCM 

19 The GEOS V2 CCM couples the GEOS-5 general circulation model (GCM) with a 

20 comprehensive stratospheric chemistry module (Bloom et al., 2005; Pawson et al., 2008). 

21 The model has 2° latitude x 2.5° longitude horizontal resolution and 72 vertical layers, with a 

22 model top at 0.01 hPa. Predicted distributions of water vapour, ozone, greenhouse gases 

23 (CO 2 , CH 4 , and N 2 0) and CFCs (CFC-1 1 and CFC-12) feedback to the radiative calculations. 

24 The performance of the GEOS V2 CCM was evaluated in detail by SPARC CCMVal (2010). 

25 The present study considers a new formulation of the GEOS V2 CCM, with an updated GCM 

26 and a new gravity wave drag scheme. 

27 

28 The present version of the GEOS-5 GCM uses the finite-volume dynamics of Lin (2004). 

29 Physical parameterizations include schemes for atmospheric convection, large-scale 

30 precipitation and cloud cover, longwave and shortwave radiation, turbulence, gravity wave 
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1 drag, as well as a land surface model. Convection is parameterized using the Relaxed 

2 Arakawa-Schubert (RAS) scheme (Moorthi and Suarez, 1992) and a scheme for the re- 

3 evaporation of falling rain (Bacmeister, 2006). RAS is a mass flux scheme with an updraft- 

4 only detraining plume cloud model and a quasi-equilibrium closure. The longwave radiative 

5 processes are described by Chou et al. (2001), and include absorption due to cloud water, 

6 water vapour, carbon dioxide, ozone, NT) and methane. The shortwave radiative scheme 

7 follows Chou and Suarez (1999), and includes absorption by water vapour, ozone, oxygen, 

8 CO2, CH4, N?0, CFC-11, CFC-12 and HCFC-22, as well as scattering by cloud water and 

9 aerosols. The turbulence parameterization is based on the Lock (2000) scheme, acting 

10 together with the Richardson-number based scheme of Louis et al. (1982). The Monin- 

11 Obukhov surface layer parameterization is described in Helfand and Schubert (1995). The 

12 Koster et al. (2000) land surface model is a catchment-based scheme that defines two soil 

1 3 layers for temperature, three soil layers for moisture, a snow pack and a canopy interception 

1 4 reservoir. 

15 

16 A previous version of the GEOS-5 GCM was used as part of MERRA, as described in 

17 Rienecker et al. (in preparation). Major changes to the moisture, turbulence and gravity wave 

18 drag schemes were made from this previous model version (Molod et al., in preparation). 

19 Increased re-evaporation of grid scale and convective precipitation, along with modifications 

20 to the turbulence parameterizations to limit the impact of the Lock scheme and enhance the 

21 impact of the Louis scheme in the presence of wind shear, improved the simulated tropical 

22 convergence and convection patterns, as well as the global stationary wave patterns. As noted 

23 by Flurwitz et al. (2010), the representation of tropospheric stationary wave patterns 

24 determines a model’s ability to simulate eddy heat flux at 100 hPa. 

25 

26 The model’s gravity wave parameterization (GWP) computes the momentum and heat 

27 deposition to the breaking of orographic and non-orographic gravity waves using the linear 

28 saturation theory by Lindzen (1981). Orographic gravity wave stress is specified using the 

29 formulation derived by McFarlane (1987) and given at the top of the subgrid-scale mountains. 

30 Subgrid-scale orography is assumed to be horizontally isotropic, and hence the orographic 

3 1 wave stress is oriented opposite to the wind direction averaged below the top of the subgrid 

32 mountains. Non-orographic wave stress is given as a Gaussian-shaped phase speed spectrum 
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1 at 400 liPa; phase speeds are assumed to be parallel to the wind direction at 400 hPa. The 

2 spectrum is composed of 9 waves with ground based phase speeds ranging from -40 to 40 m 

3 s- 1 at an interval of 10 m s' 1 , as in Garcia and Solomon (1985). As non-orographic gravity 

4 waves often accompany precipitation (e.g„ convective and frontal systems; see Richter et ah, 

5 2010), the latitudinal structure of the gravity wave spectrum is designed to mimic the 

6 structure of the climatological mean precipitation field: The spectrum maximizes (7.7><10 3 N 

7 m 2 ) at the equator, has two secondary maxima of 2 > 10 N m 2 at 60°N and 60°S, and is 

8 weakest at subtropical latitudes. The tropical peak in non-orographic gravity wave stress is 

9 necessary for the generation of an internal QBO in the model (as in Ziemke et al., 2010; see 

10 also Section 3.1). 

11 

12 In the Lindzen-based scheme, the horizontal wavelength is 100 km both for orographic and 

13 non-orographic waves. A 700 km wavelength is used for the tropical non-orographic waves 

14 to prevent an excessive downward propagation of the semi-annual oscillation into the lower 

15 stratosphere, and thus contamination of the QBO signal. The intermittency factor is doubled 

1 6 for orographic waves south of the 40°S, based on the observations of strong mountain waves 

17 from the Antarctic peninsula (Alexander and Teitelbaum, 2007) and isolated small islands 

18 (Alexander et al., 2009). This increased intermittency factor is effective in producing a 

19 reasonable evolution of the breakup of the Antarctic vortex. 

20 

21 Heat transfer due to gravity wave breaking is required for the conservation of grid-scale 

22 energy across the GWP. Heat transfer is computed from the deposition of the gravity wave 

23 energy flux into the mean flow following Warner and McIntyre (2001) and Shaw and 

24 Shepherd (2009). In addition, for the conservation of angular momentum and energy, gravity 

25 waves stress and energy flux is gradually dissipated in the top five model layers, as suggested 

26 by Shaw and Shepherd (2007). 

27 

28 2.3 GEOS V2 CCM simulations 

29 Two, 50-year time-slice simulations with the GEOS V2 CCM will be considered in Section 3. 

30 Each simulation is forced by a distinct set of SST and sea ice climatologies, each with a 

3 1 repeating annual cycle, with conditions composited from a number of observed WPEN and 
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1 ENSON events (see Table 1). Each event spans from the July preceding the SONDJF peak in 

2 tropical SST anomalies through June of the following year. HadJSSTl SSTs and sea ice 

3 concentrations at !°xl° resolution (Rayner et ah, 2003) are used to prepare the composites. 

4 The SSTs used in the WPEN and ENSON simulations are significantly different throughout 

5 most of the tropical Pacific. 

6 

7 The set of events used to create the SST and sea ice boundary conditions in the GEOS V2 

8 CCM simulations is slightly different from that used in the MERRA analysis (Table 1). 

9 HadrSSTl data were not available for the 2008 ENSON and 2009 WPEN events. The 1991 

1 0 and 1 994 WPEN events were associated with large SST anomalies in the central Pacific (Kug 

1 1 et ah, 2009) and are thus used to generate a strong model response. 

12 

13 Both simulations used fixed greenhouse gas and ozone-depleting substance boundary 

14 conditions representative of the year 2005. Variability related to the solar cycle and volcanic 

15 eruptions is not considered. 

16 

17 3 Results 

18 3.1 Climatology of the ENSON simulation 

19 The mean annual cycles of temperature and zonal wind at 60°S, in the ENSON simulation, 

20 are shown in Figure 1. The GEOS V2 CCM (Figure lc) is generally able to simulate the 

21 observed stratospheric temperature pattern (Figure la). As has been noted in evaluations of 

22 earlier formulations of this model (Pawson et al., 2008; SPARC CCMVal, 2010), polar 

23 stratospheric temperatures are biased high in the lower stratosphere in mid-winter. Similarly, 

24 the model is broadly able to simulate the SH polar jet. Modelled winds (Figure Id) remain 

25 weaker than observed (Figure lb) in the middle stratosphere, in mid- to late winter. The 

26 simulation of the breakup of the Antarctic vortex is comparable to that reported by Hurwitz et 

27 al. (2010). 

28 

29 A timeseries of zonal winds in the equatorial region, in 10 years of the ENSON simulation, is 

30 shown in Figure 2b. The MERRA timeseries of zonal winds in the 1990s is shown for 
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1 comparison (Figure 2a). The modelled QBO has a realistic period (30 ± 3 months, at 30 hPa) 

2 though the descent of easterly and westerly anomalies from the upper stratosphere is slower 

3 than observed. The modelled QBO amplitude is well simulated around 1 0 hPa but is weaker 

4 than observed in the lower stratosphere. Note that the QBO signal is absent at 50 hPa: zonal 

5 winds are generally easterly. For this reason, 30 hPa zonal winds in the 10°S-10°N region are 

6 used to characterize the phase of the modelled QBO. As for the MERRA reanalysis, QBO-E 

7 years are identified when the modelled November/December mean QBO index is less than -2 

8 m s' 1 while QBO-W years are identified when the QBO index is larger than 2 m s' 1 . For each 

9 simulation, there are approximately 25 years designated as QBO-E and 15 years designated as 

10 QBO-W. 

11 

12 3.2 Stratospheric response to WPEN events 

13 October/November eddy heat flux at 40-80°S, 100 hPa is a measure of the planetary wave 

14 energy entering the SF1 polar stratosphere (Newman et al., 2001). In the MERRA reanalysis 

15 (see Table 2), eddy heat flux in the WPEN/QBO-E composite is significantly different than 

16 that in all three of the other composites at the 95% confidence level, in a two-tailed t-test with 

17 unequal variance. Though October/November eddy heat flux in the GEOS V2 CCM (Table 

18 2) is weaker than observed, the model is able to capture some of the observed variability. For 

19 example, eddy heat flux values in the WPEN/QBO-W and ENSON/QBO-W composites are 

20 different at the 90% level. Differences between the modelled WPEN/QBO-E and 

21 WPEN/QBO-W composites are not statistically significant. As in the MERRA reanalysis, the 

22 magnitude of the eddy heat flux in the ENSON simulation does not depend on the phase of 

23 the QBO. 

24 

25 Enhanced eddy heat flux in October/November is associated with a warming of the Antarctic 

26 lower stratosphere in November/December (Flurwitz et al., 201 1). In MERRA, the Antarctic 

27 lower stratosphere warms 3-5 K in response to WPEN events in November/December (Figure 

28 3a), whiie the upper stratosphere cools by approximately 1 K. In the Antarctic lower 

29 stratosphere, the QBO-related temperature response during WPEN (Figure 3b) is nearly 

30 identical to the WPEN response itself: a 4-5 K. warming that is significantly different from 

31 zero at the 95% confidence level. That is, this warming signal can be attributed to 

32 WPEN/QBO-E events. 
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2 The GEOS V2 CCM is able to simulate the observed Antarctic temperature response to 

3 WTEN events, but not the modulation of this response by the QBO. November/December 

4 temperature differences between the WPEN and ENSON simulations, during QBO-E, are 

5 approximately 2-4 K in the lower stratosphere with a roughly 1 K cooling response above 10 

6 hPa (Figure 3c). Figure 3d shows that the modelled temperature response to WPEN is 

7 insensitive to the phase of the QBO: the lower stratosphere warms in both QBO-E and QBO- 

8 W years. 

9 

10 3.3 Tropospheric response to WPEN events 

1 1 Several tropospheric diagnostics are used to investigate the model’s insensitivity to the phase 

12 of the QBO: differences in OLR, the upper tropospheric streamfunction and calculations of 

13 Rossby wave energy in the tropical and subtropical south Pacific. As discussed by Hurwitz et 

14 al. (2011), SON seasonal mean upper tropospheric streamfunction differences illustrate the 

15 planetary wavetrains generated in response to El Nino events. Figure 4 compares the 250 hPa 

16 streamfunction in the WPEN versus ENSON cases, as well differences between QBO-E and 

17 QBO-W during WPEN events. In the MERRA reanalysis, southeastward-propagating 

18 wavetrains are seen in response to both WPEN (Figure 4a) and the QBO (Figure 4b). The 

19 model simulates the wavetrain response to WPEN (Figure 4c), with comparable magnitude 

20 and propagation direction, but lacks sensitivity to the phase of the QBO (Figure 4d). 

21 

22 Vera et al. (2004) identified a planetary wave source region in the south central Pacific, 

23 during austral spring, in WPEN-like events as compared with cold tongue El Nino and ENSO 

24 neutral events. The authors determined that the increase in planetary wave activity was the 

25 result of enhanced upper level divergence. That is, enhanced convection in the south central 

26 Pacific leads to an enhancement in Rossby wave activity and thus to a stronger poleward- 

27 propagating wavetrain in the SH. In Figure 5, OLR differences illustrate changes in 

28 convection in MERRA and in the GEOS V2 CCM simulations. Negative OLR differences 

29 imply relatively higher cloud tops and thus deeper convection. Figures 5a and 5c show the 

30 convective response to WPEN events, during QBO-E. The largest OLR differences in 

3 1 Figures 5a and 5c occur in the equatorial region, and are associated with well-established 
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1 ENSO-related longitudinal shifts in convection. The convective response to WPEN events 

2 extends to the SH subtropics: there are negative OLR differences in the Rossby wave source 

3 region identified by Vera et al. (red boxes) in both the NOAA interpolated OLR dataset 

4 (Figure 5a) and in the model simulations (Figure 5c). 

5 

6 In MERRA, there is a distinct region of negative OLR differences in the Rossby wave source 

7 region in response to the QBO, during WPEN events (Figure 5b). This result suggests that 

8 QBO-related changes in tropical circulation modulate the depth of convective activity in the 

9 subtropics, consistent with the hypothesis of Collimore et al. (2003). In contrast, there are no 

10 significant differences in OLR in this region in the model (Figure 5d). 

11 

12 Similarly, in the MERRA reanalysis and during the SON season, Rossby wave source 

13 calculations (based on Jin and Hoskins, 1995) show increases in the Vera et al. (2004) 

14 planetary wave source region (red boxes) in WPEN events relative to ENSON (Figure 6a) and 

1 5 WPEN/QBO-E relative to WPEN/QBO-W (Figure 6b). In the GEOS CCM, there is a relative 

16 Rossby wave source in the central Pacific south of 20°S during WPEN events (Figure 6c); 

17 however, this response lacks sensitivity to the phase of the QBO (Figure 6d). Possible 

18 reasons for the model’s lack of sensitivity to the QBO phase will be discussed in Section 4. 

19 

20 4 Summary and discussion 

21 WPEN is an emerging type of El Nino event (Kug et al., 2009). CCMs used to predict ozone 

22 recovery and the 21 st century stratospheric climate need to capture the extra-tropical 

23 stratospheric response to WPEN, since these events are predicted to occur more frequently in 

24 future (Yeh et al., 2009; Xie et al., 2010). 

25 

26 This study was the first to examine the modelled response of the Antarctic stratosphere to 

27 WPEN, as well as the modulation of this response by the QBO. Two time-slice simulations, 

28 one representing WPEN conditions and the other ENSON conditions, were used to both test a 

29 new formulation of the GEOS V2 CCM and confirm the observed atmospheric response to 
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1 recent WPEN events (as described by Hurwitz et al., 2011). The MERRA reanalysis was 

2 used to validate the modelled response to WPEN. 

3 

4 WPEN events enhanced poleward planetary wave activity in the central South Pacific during 

5 austral spring. In both MERRA and in the GEOS V2 CCM, enhanced eddy heat flux 

6 hastened the breakup of the Antarctic vortex by approximately 5 days at 60°S, 10 hPa and led 

7 to a relative warming of the Antarctic lower stratosphere in November/December. During 

8 QBO-E, the GEOS V2 CCM simulated a wanning of 2-4 K in the Antarctic lower 

9 stratosphere, in the WPEN simulation relative to ENSON. The modelled temperature 

1 0 response to WPEN events had the same vertical and meridional structure but was somewhat 

1 1 weaker than that in the MERRA reanalysis (4-5 K). The relatively weaker model response is 

12 likely a consequence of the experimental design: This study compared the transient 

13 atmospheric response (i.e., in the MERRA reanalysis) with time-slice simulations, each 

14 representing a near-equilibrium climate. 

15 

16 Analysis of MERRA demonstrated the sensitivity of the WPEN response to the phase of the 

17 QBO, in the 1979-2009 period. The lower stratospheric WPEN-related warming signal was 

18 as large as the QBO-related warming signal, suggesting that this signal could be attributed to 

19 WPEN events coincident with the easterly phase of the QBO. Tropospheric diagnostics 

20 suggest that it is the combination of the placement of convective activity in the subtropical 

21 Pacific during WPEN events, and the enhancement of this convective activity during QBO-E, 

22 that lead to changes in planetary wave activity in the SH and thus to the stratospheric 

23 temperature response. 

24 

25 In the GEOS V2 CCM simulations, the Antarctic temperature response to WPEN events was 

26 insensitive to the phase of the QBO. OLR, streamfunction and Rossby wave source 

27 diagnostics demonstrated that the modelled QBO did not extend far enough into the lower 

28 stratosphere and upper troposphere to modulate convection and thus planetary wave activity 

29 in the south central Pacific. In the model formulation examined in this study, zonal winds in 

30 the lower stratosphere (i.e., 50 hPa) were almost always easterly, regardless of the phase of 

3 1 the QBO at higher altitudes (i.e., 30 hPa). Thus, the model produced an easterly QBO-like 
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1 response to WPEN in both phases of the middle stratospheric QBO. If the QBO westerly 

2 anomalies were to descend further into the lower stratosphere and/or or upper troposphere, 

3 better matching observations, the model would likely be able to simulate QBO-related 

4 modulation of subtropical convection and in turn the dependence of the WPEN response on 

5 the QBO phase. 

6 
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1 Table 1. List of ENSO events composited to create the SST and sea ice boundary conditions 

2 for the WPEN and ENSON simulations. The years shown refer to the SON season included 


3 for each event. * symbols denote QBO-E years; A symbols denote QBO-W years. 
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1 Table 2. October/November mean eddy heat flux magnitude (K m s' 1 ) at 40-80°S, 100 hPa ± 

2 2 standard deviations. 
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Figure 1. Temperature (K) (a, c) and zonal wind (m s' 1 ) (b, d) at 60°S, as a function of month 
and altitude, in (a, b) a composite of ENSO neutral years in the MERRA reanalysis and (c, d) 
the ENSON simulation. 
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2 Figure 2. Equatorial (4°S-4°N) zonal mean zonal wind (m s' 1 ), as a function of time and 

3 altitude: (a) 1990 through 1999 in MERRA; (b) 10 years of the ENSON simulation. White 

4 contours denote the zero wind line. The black dotted line highlights the 50 hPa level. 
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Figure 3. November/December temperature differences (K) in (a, b) MERRA and (c, d) the 
GEOS CCM simulations. The response to WPEN events during QBO-E is shown in (a) and 
(c). The response to the QBO phase during WPEN events is shown in (b) and (d). White 
contours indicate zero difference. Black Xs indicate regions where differences are significant 
at the 95% confidence level. 
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Figure 4. Longitude-latitude contour plots showing SON seasonal mean streamfunction 
differences at 250 hPa (10' h nr’ s' 1 ) in (a, b) NCEP and (c, d) the GEOS CCM simulations. 
The response to WPEN events during QBO-E is shown in (a) and (c). The response to the 
QBO phase during WPEN events is shown in (b) and (d). White contours indicate zero 
difference. Black Xs indicate regions where differences are significant at the 90% confidence 
level in the MERRA reanalysis and 95% confidence level in the GEOS CCM. Red arrows 
indicate the approximate direction of planetary wave propagation. 
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Figure 5. Longitude-latitude contour plots showing SON seasonal mean OLR differences (W 
m" 2 ) in (a, b) NOAA and (c, d) the GEOS CCM simulations. The response to WPEN events 
during QBO-E is shown in (a) and (c). The response to the QBO phase during WPEN events 
is shown in (b) and (d). White contours indicate zero difference. Black Xs indicate regions 
where differences are significant at the 90% confidence level in the NOAA interpolated OLR 
dataset and 95% confidence level in the GEOS CCM. Red boxes indicate the approximate 
Rossby wave source region found by Vera et al. (2004). 
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2 Figure 6. Longitude-latitude contour plots showing SON seasonal mean differences in the 

3 Rossby wave source (10 10 s' 2 ), as described in the text, in (a, b) MERRA and (c, d) the GEOS 

4 CCM simulations. The response to WPEN events during QBO-E is shown in (a) and (c). The 

5 response to the QBO phase during WPEN events is shown in (b) and (d). White contours 

6 indicate zero difference. Black Xs indicate regions where differences are significant at the 

7 95% confidence level. Red boxes indicate the approximate Rossby wave source region found 

8 by Vera et al. (2004). 


